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A qubit (a spin-1/2 particle) prepared in the up state is scattered by local spin-flipping potentials 
produced by the two target qubits (two fixed spins), both prepared in the down state, to generate an 
entangled state in the latter when the former is found in the down state after scattering. The scat- 
tering process is analyzed in three dimensions, both to lowest order and in full order in perturbation, 
■ with an appropriate renormalization for the latter. The entanglement is evaluated in terms of the 

' concurrence as a function of the incident and scattering angles, the size of the incident wave packet, 

, and the detector resolution, to clarify the key elements for obtaining an entanglement with high 

^Sj ■ quality. The characteristics of the results are also discussed in the context of (in)distinguishability 

of alternative paths for a quantum particle. 

PACS numbers: 03.67.Bg, 03.65.Nk, ll.lO.Gh 

I. INTRODUCTION 

S'- n . . . . 

O • Entanglement plays a crucial role in the field of quantum information and technology [![, though its acquisition or 

controlled generation is by no means a trivial matter: it is one of the most peculiar features of quantum theory with no 
I classical analog, and there are several proposals for its generation. One often makes use of their mutual interaction to 
d ' make the two quantum systems entangled On the other hand, when they are separated far away and/or when their 
^ , mutual interaction is considered absent or negligibly weak, one may resort to a third quantum system to make them 
^ , ' entangled through its individual interaction with each of them d, [j, H, H, 0, 9 ■ For example, one can consider single 

modes in two cavities with no direct interaction two spatially separated atoms 0], two remote atomic qubits 
OA • or a system of two magnetic impurities (spins) embedded in a solid 0, ^'^d regard an atom, a single cavity mode, 
^ " a photon, and an electron spin, respectively, as mediators of interaction between distant quantum parties, in the 
laboratories. This kind of scheme with the use of "entanglement mediators" has been investigated for simple systems 
of qubits (quantum two- level systems) , usually with the assumption that the strength of the interaction between the 
mediator and each quantum systern, i.e., the magnitude of the coupling constants and the interaction durations, are 
completely under our control 0, |3j ■ Conditions under which maximally entangled states are realized are expressed 
in terms of coupling constants and interaction durations in these cases. 

Even though the assumption of complete controllability of such experimental parameters as interaction strength 
[ is considered to be legitimate, for example, when the interaction is well controlled by switching on/off the external 
parameters [H, there are still cases in which such an assumption is untenable or its applicability is questionable. In 
particular, when the interaction time is not well defined or its definition necessarily requires a resolution in some 
conceptual issues, like the definition of moments of the beginning and the end of interaction for a particle described 
^ ' by a wave packet with a finite width and scattered by a static potential, we would be forced to treat the process as 
5^ , a quantum mechanical scattering process of a mediator system off the target, where additional (internal, e.g., spin) 
degrees of freedom are duly taken into account. In scattering processes, such notions as the first or the last moment 
of interaction are not considered to be proper issues to be asked, for the interaction is supposed to be gradually 
turned on and off and the scattering matrix, which describes the transition from the remote past {t —^ — oo) to the 
remote future [t — > -l-oo), could be the only quantity of physical relevance. The interaction strength is in a sense 
automatically and implicitly given and we have no choice of defining or controlling the interaction duration once the 
initial conditions have been fixed. It is therefore an interesting and nontrivial matter of physical relevance to examine 
whether the schemes for entanglement generation or extraction, based on the interaction between the mediator and 
subsystems, could remain valid and function properly even when one has little controllability on such parameters as 
time. 

In this paper, a three-dimensional scattering process, in which a qubit (a spin-1/2 particle, playing a role of 
mediator) is scattered off a fixed target composed of two qubits by spin-fiipping delta-shaped potentials, is considered 
to examine the ability of obtaining an entangled state in the target system when a spin flip has been confirmed in the 
final state of the mediator qubit. The same setup has already been considered, but essentially only in one (spatial) 
dimension, to generate a maximally entangled state in the two-qubit system 0, 13 • It is shown that a maximally 
entangled state is obtained if one can properly tune the interaction strengths. On the other hand, if it is treated 
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as a scattering process in one (spatial) dimension, the entanglement can be enhanced by a resonant scattering when 
the incident momentum of the mediator (or the distance between two target qubits) has been properly chosen in the 
initial setup [1]. Notice that the treatments of the scattering processes so far are not considered to be completely 
satisfactory because in one dimension, there would be no way to incorporate such important physical parameters as 
incident and scattering angles, the lateral size of the wave packet and the detector resolution (e.g., an aperture of 
detector mouth). The purpose of this paper is to take these elements into account in the scattering process and to 
clarify the dependence of the resulting entanglement on these parameters. 

The paper is organized as follows. The Hamiltonian for our three-dimensional scattering process is presented in 
section Sec. |lTl The scattering matrix (S matrix), which describes the transition of the mediator qubit from a given 
initial state to the final state, is then introduced in Sec. lIIII If the scattered qubit is detected with its spin flipped, one 
knows that the target system of two qubits is in an entangled state, provided that the both target qubits have been 
prepared in the down state and the interaction preserves the total number of spins-up. The degree of entanglement is 
measured in terms of the concurrence and the concurrence is expressed in terms of the S-matrix elements. In Sec. IIVI 
the S-matrix elements and the concurrence are calculated to the lowest order and the dependence of the entanglement 
on physical parameters and its characteristics are discussed. In this case, the (in)distinguishability of mediator's 
alternative paths to the detector is shown to be closely connected with the values of concurrence and therefore with 
the degree of entanglement obtained. When one tries to evaluate higher-order contributions and goes to higher- 
order terms in perturbation theory, one realizes that the delta-shaped potentials in three dimensions bring about 
(ultraviolet) divergences and a proper treatment of such divergences (renormalization) is required @, We follow 
the prescription proposed by Jackiw [l^l to treat delta-shaped potentials and finite S-matrix elements are obtained 
in Sec. |V] in full order in perturbation theory by properly introducing counter terms. The resulting concurrence is 
examined and its characteristics are discussed in comparison with its lowest-order counterparts. The final Section IVll 
is devoted to the summary and outlook. The details of the calculations that are too involved to be presented in the 
text and some related aspects arc shown in Appendices |X] and |B] 

II. SETUP OF THE PROBLEM 

The setup is sketched in Fig. [T] We try to make two qubits A and B, both prepared in the down state HDab, 
entangled. Notice that the initial state of A+B is separable and not entangled at all, HDab = |i)A ^ IDb- Assume 
that these qubits are fixed at positions — d/2 and d/2, and there is no mutual interaction between them. In order to 
make them entangled, another qubit X is prepared in the up state |'f)x and is projected on qubits A and B. Qubit 
X is then scattered by spin-flipping delta-shaped potentials produced by qubits A and B and is finally detected by a 
spin-sensitive detector. This physical process can be described by the total Hamiltonian 

i7 = Pl+g53(r + d/2)(^(^) .^(A)) +g<53(r-d/2)(^W •^(B)N (3.!) 
2m 

where p and r are the momentum and the position of qubit X of mass m and cr^''' (J = X, A, B) the Pauli matrices 
acting on qubit J. The interaction between X and A(B) is described by the spin-exchange operator multiplied by 
a short-range (delta-shaped) potential with strength g and the energy difference between the up and down states is 
assumed to be neglected (or the same for all qubits and therefore neglected). We shall treat two qubits A and B 
symmetrically, for simplicity. 




FIG. 1: (Color online) Qubit X is scattered by qubits A and B and is detected by a detector with a finite resolution. 
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Since the Hamiltonian H preserves the total number of ups (t) among the three qubits X, A, and B. if we find qubit 
X in the spin-flipped (i.e., down) state at the detector after scattering, we are sure that one of the qubits, A or B, 
must be in a spin-flipped state, that is, system A-l-B is either in |ti)AB or in |it)AB- If there is no way to judge which 
spin has been flipped during the scattering by qubit X, the state of qubits A and B is certainly in their superposed 
state, that is, an entangled state, like a|ti)AB + /3|iT)AB- This argument is valid when the interaction strength can 
be freely adjusted by some external parameters. Indeed, it has been shown that a maximally entangled state can be 
extracted by properly tuning the interaction strength, that is, [the coupling constant] x [interaction duration (time)]. 

It is, however, to be noticed that when the spatial degrees of freedom of the particles and local potentials are 
incorporated, as in (|2.ip . to describe the interactions between qubits X and A(B) as a scattering process, the situation 
would change drastically: The state finally be extracted can no longer be a pure state, because it shall be given 
as a reduced density matrix after being traced over possible spatial or momentum degrees of freedom of X within 
the detector resolution and therefore it becomes mixed in general. Furthermore, it is not at all trivial whether we 
can still expect an entangled state with high quality to be extracted, since we would have no direct control on the 
interaction time in the scattering process. (We assume adiabatic switchings of interaction and no (abrupt) changes 
of the parameters in the Hamiltonian are considered to occur.) Even if we could deliberately choose the incident 
momentum of qubit X, so as to adjust effectively the interaction strength, since a quantum particle with a less 
momentum spread is necessarily accompanied with a longer wave packet in space and the interaction region or range 
cannot always be given precisely (except for some ideal cases) , it is almost impossible or meaningless to talk about the 
precise moments of particular events in the scattering process. Stated differently, the quantity of physical relevance 
in the scattering problems is the scattering (S) matrix elements, which describe the transition from the remote past 
t —oo to the remote future t +00 under the Hamiltonian H. This is one of the motivations of the previous 
works where the ability of extraction of entanglement has been examined in the context of the quantum mechanical 
scattering process in one spatial dimension [3, § . Here the problem is generalized and extended to three spatial 
dimensions, where a much more rich variety of physical parameters are expected to play roles. 

Let the incident qubit X be described by a Gaussian wave packet, with a central momentum hka and a spatial 
width w. Qubit X is sent exactly towards the midpoint of target qubits A and B, and the normalized initial state of 
the total system reads as 

// O 2 \ 3/4 

where |fcTii)xAB = j^Dx <8) |ii)AB- A spin-sensitive detector, placed far away from the target detects qubit X 
scattered in the direction Aid seen from the origin (scattering center) with an opening angle A9. The detector is 
assumed to be indifferent to other quantities than spin degrees of freedom of qubit X and we are interested in only 
those events in which scattered qubit X has been found in the spin-flipped (i.e., down) state |i)x- When qubit X has 
been found in ||)x in the detector that covers the solid angle 

Ail^ dip deisin^ = 47rsin2(Ae'/2) (2.3) 

around direction fco, the state of the target system A-)-B is given by the reduced density matrix of the form 

/O 0\ 

P=^l^^d^kM\S\^o){MS^\HU = ^l 21 , (2.4) 

\0 0/ 

where S denotes the S matrix and we take |tT)AB! |Ti)AB! |iT)AB! and |ii)AB as the standard basis for the matrix 
representation of the density operator p. Observe that owing to the conserved quantity, i.e., total number of ups, only 
those elements relevant to |ti)AB and |J,t)AB can have nonvanishing values in the current setup. The normalization 
constant 

P = tr(/ rf3fcx(fci|^|^o)(V'o|^^|fci)x| =aii+a22 (2.5) 

I J AO J 

is nothing but the probability that this particular event occurs, i.e., the yield. Given the state p for two qubits, the 
degree of entanglement can be measured in terms of its concurrence C{p), which reads, for the above p in (|2.4[) . as [ll| 



(2.6) 
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III. SCATTERING MATRIX 



We proceed to the calculation of the matrix elements a-y in p.4p . The scattering matrix S is defined in the following 
limits 

lim e*^°*/'^e-'"^(*-*')/''e-'■^°*'/^ (3.1) 

t — !- + 00 

t' — > — oo 

where Hq = /2m is the free Hamiltonian. Let |fcC) be an eigenstate of iJo with C collectively denoting the spin 
state of three qubits X, A, and B, i.e., 

H^\kO=Ek\k(:), Ek = ^, (3.2) 

Zm 

and I'I'fcC) an eigenstate of the total Hamiltonian H = Hq + V belonging to the same energy Ek- Notice that since 
the free Hamiltonian Hq is independent of the spin degrees of freedom C, the cigenstates are degenerated with respect 
to C. The normalized solution of the eigenvalue equation 

H\^kO ^Ekl^kC) (3.3) 

is formally given by 

1^.0 = IK) + ^;-^yH) 

(*fcCI*fc'C'> ='53(fc-fe')<5cc'. (3.5) 

Its coordinate representation reads as 

{r\^kO = {r\kO - I d^r' Gk{r - r' )'^V {r' ) {r' \^ , (3.6) 
where the retarded Green function Gk is given by 

^"^"■^ ^2^ J (27r)3 E,~Ek-ie^ i^'' ^^'^^ 

(The above eigenstate satisfies the so-called "out-going wave" boundary condition for r oo.) These eigenstates 
{l^fcC)} form, together with possible bound states {|?ix)} with negative discrete energies Enx < 0, a complete set 

I d'k |vI/fcC>(*fcCI + E \^x){nx\ = 1- (3.8) 

C n,x 

We are now ready to evaluate the S-matrix elements. Notice first that the following limits are evaluated by making 
use of the above eigenstates l^'fcC)? 



C 



/ '^'fc'|vl/fe,C')(*fc'C'IK)e*<^''-'"^'='*/''+5]|nx)(nx|fcC)e 
*fcC) as t ^ -oo, 

*feC)+27rz^ /"d3fc'J(£;fc-£;fc,)|*fc'C')(*fc'C'I^^IK) as t^+oo. ^^'^^ 

These results correctly reflect the fact that the eigenstate \^kO satisfies the "out-going wave" boundary condition 
and no bound state can survive in the asymptotic regions where the energy conservation is recovered. The S-matrix 
element is thus given by 

{k'C\S\kO = S^k' - fc)<5c'c - 2Tnd{Ek' ~ -Bfc)(fc'C'|F|*fcC)- (3-10) 
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IV. BORN APPROXIMATION 



To lowest order in the coupling constant g, the eigenstate l^'fcC) in the S-matrix element (|3.10[) can be replaced 
with its free counterpart \kQ since V is proportional to g. Up to the first order in (7, the relevant S-matrix elements 
read as (subscripts x,a,b for the spin states shall be suppressed for notational simplicity in what follows, provided no 
confusion would arise) 



{km\SH'o)=9Aik) + 0{g'), 
{km\S\i'o) = -gA*{k) + 0{g') 



(4.1) 



where 

Aik) = J <fk' Mk')S{Ek - Ek')e^''-'''^-'"\ (4.2) 

The matrix elements aij are simply expressed as 

an = a22 = f d^k \A{k)f ee a, 012 = -g^ [ (fk A^{k) ee a, (4.3) 
JAn JAn 

and the concurrence C and the yield P are given by 

C{p) = p ^ 2a. (4.4) 

a 

When the incident wave packet ipoif^) is well monochromatized, wko 3> 1, the above quantity A(k) is approximately 
evaluated analytically. Since it is expected that the main contributions in the amplitude A{k) in (|4.2p . 

A{k) = —('^^^^ ^ 1 ^Z^^-w^-q^ ^l{k-kQ-q)-d/2^is{Ek-Eka-ti^ka-q/m~E,) 

27r2 V TT y 27r i 

are due to the integrand with small \q\ < l/w, we would be allowed to drop the term Eq in the third exponent relative 
to EkQ, for Eq < /2raw'^ <C h^h:^/2m = Ek^- Then, the integrations over q and s are easily performed and we arrive 
at 



The matrix element a = an = 022 is now reduced to 

3/2 



a^^(^) ^]^,-ld^-ikrdf]/s..-^ f ^^^^-^^uUEUE.,-if/2, (4.7) 



The integration over k can be estimated as 



An 



2m^El 



00 



^E,^ eVB"-^^-^ . An ^E,,, (4.8) 



in the lowest order in 1/w k^. Thus, we obtain 



We apply the same approximation to the other matrix element 5, = ai2 to get first 
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a • 



.g2 /27r 



3/2 



mw 



47r4 V J 27rh?Ekr, 



Ko JAQ 

The last integration is similarly approximated as 



-[d-'-{ko-df]/Sw^ / ^3f^^~w^kl{Ek/Ek„-lf/2^-iko-d{Ek/Ek„-l)/2^i{k-ko)-d_ (4.10) 



^3j^^-w^k^iEJEkg-lf/2^-iko-d{Ek/Ekg-l)/2^Hk-ko)-d 



AO 



(fit / (;;2,g-«'"'=52;^/2-i(fco-d)2:/2+4fco(fc-d)/T+^+(l/2)ln(l+2;)g-ifco-d 
Af2 J-1 



Ej,.„ I ffl.p-[{kok-koyd]''/Sw^kl^i(k„k-kayd 

An 



which leads to 



q - , ^-\d''-(ko-df]/Sw' f ^2^^-[ikok-ko)-df/8w'k'a^i{kak-ko)-d 

We finally end up with the following expressions for the concurrence and the yield 

lal 1 



a An 



^2^ ^ikok-d-[{k^ko)-df /8w^ 



An 



(4.11) 



(4.12) 



(4.13) 



P = 2a- 



2 2 



(4.14) 



The angle integrations over fc in the concurrence C are numerically performed and the result may be plotted as a 
surface in three dimensions, swept by a vector, originating from the scattering center and pointing to the detector, 
i.e., representing the scattering direction of qubit X or the detector direction fc^, with its magnitude equal to C. 
The following figures show sections of such surfaces, which are cylindrically symmetric with respect to vector d, for 
different values of parameters. 



A. Characteristics 



1. Cases when the incident wave packet is almost like a plane wave: w ^ d 

Figure [2] shows the concurrence C generated by the incident wave packet with a large spatial width w ^ d. In 
this case, (i) the concurrence C is independent of the incident direction ■ d = cosOq. It is clear from these figures 
that (ii) it depends on in which direction the scattered qubit X is detected (dependence on the scattering direction 
fcD • d, ~ cosOd) and it takes the maximal value when qubit X is captured in the direction of id, i.e., on the line 
connecting target qubits A and B. The opening angle A9 of the detector mouth (the detecting resolution) also affects 
the concurrence C: we understand that (iii) the concurrence C is reduced considerably when A9 is increased, while it 
still keeps the maximal value in the id directions. Finally, (iv) the concurrence C is small for a large incident wave 
momentum fco or a short wavelength compared with the distance d between the two qubits A and B in the target. 



2. Cases when the incident wave packet is small: w < d 

If the incident wave packet is small compared with the distance between two qubits A and B in the target, i.e., 
w < d, (v) the concurrence C becomes smaller in general and (vi) the dependence on the incident angle ^o, which is 
almost absent for w 3> in Fig. [21 appears as seen in Fig. [31 We observe however that (vii) there are directions where 
the concurrence C takes the same value as that for the case with a large w shown in Fig. [21 that is, when qubit X is 
detected in the same direction as the incident direction fcp or ii^ its cylindrically symmetric directions with respect to 
d. 
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FIG. 2: (Color online) Concurrence C in the Born approximation for cases with a large incident wave packet, w ^ d, shown 
in polar coordinates with radius C as a function of the scattering angle 9u relative to the alignment d of the target qubits A 
and B. Parameters are: (a) A6 — it/25, n/15, n/6 with kod = 10 (dependence on the opening angle of the detector mouth 
A9), and (b) kod = 5, 10, 25 with A6 = n/15 (dependence on the incident momentum fco). The thin-line circles indicate C — 1 
(maximal entanglement). 




FIG. 3: (Color online) Concurrence C when the incident wave packet is small, w/d = 0.5, is shown for different incident angles 
6*0 = 0°, 60°, 90° relative to the alignment d of qubits A and B, with kod = 10 and A6 — 7r/15. The concurrence C for a 
large incident wave packet w/d = 10 (with the other parameters the same) is also shown (dashed curve) as a reference. In the 
directions indicated by the dots (where 6u = Oo, i.e., in the same direction as the incident direction fco and in its cylindrically 
symmetric directions with respect to d), the concurrence C takes the same values in both cases ui ^ d and w < d. 



3. Yield P 



According to the mathematical expression (|4.14p . the yield P is independent of the scattering angle O-q: no Young- 
type interference pattern is observed under the condition of the detection of X with its spin flipped. On the other hand, 
it depends on the incident direction feg • d and becomes maximal when qubit X is injected along the line connecting 
A and B, i.e., fco || d. This may be due to the effective density of the incident probability current felt by the target 
qubits and the dependence on the incident angle disappears when the incident wave packet becomes broad enough 
w ^ d so that it is seen from the target uniformly spread in space. 
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B. Conditions for Obtaining Higher Entanglement 



The above characteristics (i)-(vii) of the concurrence C can be understood on the basis of its mathematical expres- 
sion (I4.13|l . For the case of large ^ c?. it could be further reduced to 



C{P) 



1 



An 



(4.15) 



In order for this quantity to be nonvanishing, the phase of the integrand fcofc • d should not change considerably within 
the integration domain Af2. That is, the condition under which higher concurrence and therefore higher entanglement 
is attainable reads as 



k^k ■ d 



3-Ae 



D+Ae 



= 2fco(isin 9^) sin A6' < 2tt for w d. 



(4.16) 



We understand that this condition well describes the characteristics (i)-(iv) mentioned above in Sec. IIV A II Interest- 
ingly, it is possible to perform the angle integrations in (|4.15p analytically (Appendix |^ , from which one can derive 
an approximate expression for the concurrence 



C{p) ~ 1 - -{kod)^ sin^ Ou sm^{A9/2) 



(4.17) 



when the opening of the detector is sufficiently small A9 <C 1. It would be evident that the above characteristics are 
again well explained by this approximate expression. On the other hand, if the incident wave packet is small u; < d, 
we have to keep the second term —[{k — ko) ■ rf]^/8w^ in the exponent in (|4.13p . Since this term would entail an 
exponential reduction of C, the condition for keeping a higher concurrence becomes 



\{k — ko)-d\^w for w < d, 
which explains well the characteristics (v)-(vii) observed in Sec. IIV A^ 



(4.18) 



C. Indistinguishability 

It would be interesting to interpret the above conditions (|4.16p and (|4.18p for higher concurrence C in the context of 
the (in)distinguishability of the paths taken by particle X. As a general rule in quantum theory, (in)distinguishability 
of alternatives results in (non)vanishing of quantum interference [T2| . Since the concurrence C is proportional to the 
absolute value of the off-diagonal matrix clement d = ai2, its value is rather dependent on the information about 
which qubit A or B has changed its spin state in the scattering process. If one could obtain such information in 
principle, there remains no quantum interference between the two alternatives and the off-diagonal elements of the 
density matrix become vanishingly small. On the contrary, if there is no way to obtain such information, one can 
expect a high quantum interference retained and hence a high entanglement in A and B. 

Notice also that in the lowest-order perturbation qubit X with its spin flipped has been scattered only and surely 
once, either by qubit A or B, and the interaction has certainly changed their spin states: there is a direct link between 
the information about which qubit has scattered X and that about which spin has been flipped. Thus we expect that 
if one could distinguish the two alternative paths of qubit X, originating from qubit A or B. one is able to know which 
spin has been flipped. This knowledge would result in a reduction of the off-diagonal matrix elements and therefore 
of the concurrence C. 

We understand that the conditions for higher concurrence (entanglement) (|4.16p and (|4.18p can be interpreted as 
those for the indistinguishability of the two alternative paths from qubit A or B. Indeed, if the incident wave packet 
is long w ^ d and therefore is approximately considered as a plane wave of wavelength Aq = 27r//so, the condition 
for higher concurrence (|4.16p is understandable in terms of the resolving power of an optical device. It is known in 
classical optics that the optical device that has an aperture A9 seen from an object composed of two optical sources 
with mutual distance d is unable to distinguish the two sources if the condition 

< (4.19) 
d sm t^D 

is satisfied (for A9 ^ 1) [l3|. See Fig. HI This is essentially the same condition as ()4.16p . On the other hand, in the 
opposite case with it; < d, since the quantity |(fc — fep) • <i| is nothing but the difference in length between the two 
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FIG. 4: (Color online) Resolving power of an optical device and the condition (|4.16p . The optical device with an aperture 
can resolve the separation d if the wavelength Ao = 27r/fco is short enough, i.e., when Ao < d AS sin 6^0. Otherwise, it cannot 
distinguish two sources A and B, the condition of which is nothing but (|4.16p . 




ikQ - d 3 —k ■ d 



FIG. 5: (Color online) Two paths via A and B have difference in their lengths by |(fco — k) ■ d\. 

paths via qubit A or B (sec Fig.[5|); if the condition (|4.18p is not satisfied, one can determine the path particle X has 
passed through on its way to the detector, for the difference in the path length is certainly larger than the size of 
the particle, w. In this way, the conditions for higher entanglement, in both cases, are interpreted as those for the 
indistinguishability of the paths taken by X. 

The one-to-one correspondence between the path of X and the flipped spin A or B also explains why no interference 
is observed in the yield P (Sec. IIV A3|) . One can know the path of X, either via A or B, by looking at the spin state 
of A and B after scattering. This accessibility to the information on the path of X erases the interference between the 
two alternative paths. 

Notice, however, that these relations of the concurrence C and the yield P with the indistinguishability are limited 
to the lowest-order results and further study is necessary for higher-order terms in perturbation theory, where multiple 
scatterings, which invalidate the one-to-one correspondence between the knowledge of the particle paths and that of 
spin flips, are present and a resonant scattering is expected to play a role. 

V. ESTIMATION OF HIGHER-ORDER TERMS 
A. Spin-flipping delta-shaped potential and renormalization 

The multiple scatterings are absent in the lowest-order terms in perturbation and their effects, including those of 
resonant scatterings, can be seen only when one goes to the higher-order calculations. In this case, one realizes that 
a proper treatment of them is required, for otherwise the result would become trivial, that is, no scattering would 
occur by the delta-shaped potentials in dimensions greater than one |9|. This is a famous issue in quantum theory and 
Jackiw has proposed a prescription how to deal with such systems jlO| . According to Jackiw's prescription, we have 
to renormalize the strength of the delta-shaped potential so that the source term becomes nonvanishing. In other 
words, the coupling constant in the Hamiltonian, g in our case, has to be regarded as a bare one and absorb possible 
divergences arising from the Green function at the origin Gk{0). 

In our case, however, another element that was absent at that time, i.e., the spin degrees of freedom, come into 
play, which would require another care in dealing with higher-order terms. In this respect, it is important to notice 
that the interaction of the form gcr^^^ ■ cr^^^ inevitably causes another type of interaction in its higher-order terms. 
For example, in its second order, a term proportional to the unit operator in spin space, that is not proportional to 
the original form of the interaction, appears: 



(5.1) 
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This is easily understood from the fact that the interaction can be written as 

^(X) . o-(A) = 7?3XA - ST'iXA, (5.2) 

where VsxA and VixA are the projection operators on the spin-triplet and singlet spaces, respectively, 

^3XA = J(3 + <tW-<t(A)), 1(i_^(x).^(A))^ (5 3) 

^3XA=^3XA, 'PiXA ^IXA, 'PsXA^'lXA = T'lXA^'aXA = 0, PsXA + PlXA = 1 • (5.4) 

This means that in any higher-order terms in pertm'bation, there are only two types of interactions proportional to 
"PsxA or VixA present and we need to renormalize these two terms simultaneously. That is, we need two counter 
terms to obtain sensible results. Let us therefore introduce another bare coupling constant g' and start from a bare 
interaction Hamiltonian 

V{r) = S'^ir + (i/2)(5(TW . cr^^) + g') + ^3^^ _ d/2)(gcr(x) . cr'^) -t- 5'), (5.5) 
which is also written as 

V{r) = S^r + d/2)[{g + g')VzxA + (-3.g + .g')^ixA] + 5\r - d/2)[{g + g')VzxB + {-ig + <?')^ixb] 

= ^(53(r + d/2)QxA + ^5\r - d/2)QxB. (5.6) 
We can formally solve the equation for the eigenstate given in (|3.6p in the coordinate representation {r\^kCl^ 

(r-l^-fcC) = {r\kC) - Gk{r + d/2)QxA(-d/2|*fcC> - Gk{r - d/2)QxB(d/2|*fcC), (5.7) 
from which the source terms oc (±<i/2|5'fcC) can be derived by solving 

/l + Gfe(0)QxB Gk{d)QxA \ / (d/2|ffcC) \ ( {d/2\kO \ 

V Gk{-d)QxB l + GkWQxAj V(-d/2|*fcC)/ " \{-d/2\ko) ' ^^'^^ 
After an elementary but a little bit tedious calculation, we formally arrive at 

QxA(-d/2|vI/fcO = ■^ i_G,(-l)BGk(d)A (^"'^/^IfcC) - Gk{-d)B{d/2\kC)) , (5.9a) 

QxB{d/2\^U) =B ^_ ((d/2|fcC) -G,(d)^(-d/2|K>), (5.9b) 

where we have defined 

It would be evident that these expressions (j5.9ap and (j5.9b[) clearly represent multiple scattering processes in terms of 
the effective (self-)couplings at qubit A and B, A and B, and the amplitudes Gk{d) and Gk{—d), describing particle's 
propagations. A— >B and B— >A, respectively. 

At this point, we have to recall that in dimensions greater than or equal to two, the Green function at the origin 
Gfc(O) is divergent. Actually in three dimensions, it diverges linearly with a cutoff A, 



d^q 1 



{2nf g2 



I , f 1 1 



lim / dqq 



A— >oo (27r)2 Jq \q — k — ie q + k + ie 

= lim -1-A + Ofe, Ofc = ^. (5.11) 

Therefore we need appropriate renormalizations of the coupling constants in order to obtain nontrivial scattering 
amplitudes, even though we also know that, according to the dimensional regularization, there are no ultraviolet 
divergences in odd dimensions and it gives exactly the same finite term. Since the above operators A and B are 
explicitly evaluated, e.g.. 
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^ = T? [(.9 + 9')'P3XA + (-3.9 + g')VixA] - ^ 



2?7i .9 + .g' -T, , -"^9 + 9' ^ 

-^3XA + -7——^ 7 — ; i^-y ' IXA 



;i2 Vl + 2?7i(5 + 5')Gfc(0)/;i2' •'-^ ■ i + 2m(-3.9 + .g')Gfe(0)/;i2' ^--y (5-12) 

and the divergences only appear through these operators A and B in (jS.lOp . the following renormalizations of the 
bare coupling constants are sufficient to make everything finite and nontrivial, 

1 1 277;. 1 . ^ 

lim T-^A, (5.13a) 



9r + 9'r 9 + 9' ^ A^oo 27r2 

1 1 2m 1 
T = 7 + TT lim TT^A, (5.13b) 

by which the above A{B) reads as 

-^(^^ = ¥ (i + 2m(t+t;)f^./ft^^'^^<^' + TTM^S4)fV^^'''^^^0 ■ ^^-''^ 



B. Multiple scatterings in terms of projection operators 

Assume that the renormalized interaction takes the spin-exchange form (|2.H) with the renormalized coupling con- 
stant gr, in other words, we shall set the other coupling constant vanishing g'. = 0. In this case, we have 

(l + W3)(l-6)^ ^ ^ 

where = drngrflk/fi^, and therefore 

a(-.)8a(.).4 = Giw (^)' (i^^,^i_^,). (-^"" ■ - - - &) 

s jHi„e^> ^ (<7l=l X <7l»l) - &<7l=" ■ (,t'»1 + <7l=l) + <7l»l ■ ,t"» + (I] (5,16a) 

and similarly 

Gkid)AGki-d)B = A^[-7^(^) • (<t(^^' X <t(A)) - ^T*^) • (^(^' + <t(s') + <t(^) . ^(^) + 4^] (5.16b) 

with 

In order to evaluate the higher-order terms, which have formally been summed up like (|5.9p . one needs to know 
the powers of the above quantities G]:{d)BA and G^{d)AB. At this point, observe that the three spin operators 
o-(X) • (cr(^) X cr^^)), (T^-^^ • (cr(^) + cr'^'^)) and cr^*^) • cr^^) form a closed algebra with respect to the symmetrized 
multiplications among them. This implies a possibility that there are projection operators made of these three 
operators, in terms of which the operators BA and AB can be expanded uniquely. Indeed, one can confirm, after 
elementary but tedious calculations, that operators defined by 



■Psil) ^ ^ J(i + (\ 7)-(^) • (-(^) X <t(A)) + U\-i) crO^^ ■ (a(^) + <t(b)) + ^a^^) . ^(b) + (5.18a) 



(5.18b) 

7^± ^ ±i[<TW • (^(A) + <t(b)) + ^(A) . ^(B)] + 1 ^ (5.18c) 
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with = 1 and 7 an arbitrary parameter, satisfy the following projective relations 

Vs{i) + Qsh)^l, Vsh)Qs{l)^Qs{l)V,{^) = Q, Vl{^) = Vs{i), (5.19) 

n+ + n^ = i, n+n^ = n^n+^o, n\ = n+. (5.20) 

Notice that the operators TZ± are not necessarily orthogonal to Vsij) and Qs{j)- In fact, 

rsh)n+ = n+Vs{i)^n+, (5.21) 

so that TZ^ is orthogonal to Qsij), but 

Qsh)n- = n^Qsh) = Qsh), (5.22) 

which means that TZ_ is not orthogonal to Vsij) and 

7',(7)7^- = n^Vsh) = Vsh) -n+ = q^, (5.23) 

that is orthogonal to Qsi'^), as it should be. The meanings of these projection operators are exposed in Appendix IB] 
Three operators Vsil), Qs{l), and TZ+ may be used to uniquely expand the relevant operators BA and AB. This 
is indeed possible and one finds that 

Gl{d)BA = aV^ij) + pn++SQ^i^), Glid)AB ^ aV±i^) + f3n+ + 5Q±{^) (5.24) 

with appropriately chosen parameters 

7=^fl±2z , ^ + ^^ ] (5.25) 

and 



« = A'(lT*V3-(l+a-)') , /? = 2/1(1 -2efc±zv/3- (1 + 6)2), <5 = /2(l±V3-(l+efe)2) • (5.26) 
The projective properties of the operators Vsij), Qs(7), and TZ+, shown in (|5.19p - (|5.2ip . easily lead us to 



and 



^ -^±(7)+^^ 7^-3T-T^V+ + T^2±(7)- (5.27b) 



1 - Gkid)AGki-d)B 1-a Vl-(" + /3) 1 - Q'/ 1 - 

C. S-matrix elements and the concurrence C in full order 

The relevant S-matrix element is now expressed as (Cab — Ti oi |t) 

(fciCABl^l^o) = I d^k' {kiCAB\S\k'ni)Mk') 

= f d^k' Mk'){-'27ri)d{Ek- Ek'){ki<:AB\V\-i'k'tii) 
d^k'ijo{k'){~27n)S{Ek~Ek') 

X ^ ((fc|-d/2)(iCAB|QxA(-d/2|vI'fc,tii) + (fc|d/2)aCAB|QxB(d/2|vl/fe,tU) 
/ ^Vo(fc')'5(£;.-i?fe') 

(^2 



- (iCAB||;^^^G.(.)^-^i^|Tii)e^(-'=>^ 
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- (iCAB||^^-4G.(.)^-^l^|Tii)e-(-'=')-/^) . (5.28) 

By plugging the explicit forms of the inverse operators in terms of the projective ones in (j5.27p and evaluating the 
matrix elements of the three spin operators, cr'-^' • (cr'^) x cr^^)), cr^^^ ■ (cr^^^ + cr''^^), and cr'^^ • cr(^\ the relevant 
matrix elements are calculated, under the same approximation as we have taken in deriving (|4.6|) . to be 

(fcin|5|^o> ~ Mu[A{k)[l + fl{l - il)] + A*(-fc)[-(l + a) + (1 - + 'i)fl]fk 

- A*{k){l - a)h(l + 6) + (1 - ikfiik + i)fl]fl - A{-k){l + - (Wk) 

= AAfe[A(fc)afc + A*(fc)e*'*fefc][l - e-*'*(l (5.29a) 

and 

(fciitl^l^o) ~ A4.{A(fe)(l - a-)[-(l + a-) + (1 - £.k?{ik + 3)/|]/| + A*(-fc)(l - a-)[l + /l(l - ek)\h 
A*{k)[\ + /,2(1 - a)] A{-k)[-{\ + Cfe) + (1 - ikfiik + 3)/2]a} 
= -AAfc[A*(fc)afc+A(fc)e-*%.][l-e*'*(l-a)/fc], (5.29b) 
where A{k) is given in (14. 6p . 



A4 



1 ft^ 



(l-a)(l-(5)[l -(« + /?)] 2mGfc(d) 



1 f , 



k 



(5.30) 



[1 - 2fk + (1 - &)(3 + + 2A. + (1 - e;.)(3 + - A'(l - 6)2] 2m Gfe(d) ' 

and 

a,. = 1 + A2(l - e^), 6fe = [-(1 + a-) + (1 - aO'(a + 3)/|]A.. (5.31) 
Here we have made use of the relation 

A{-k) = e-'''-'^A{k). (5.32) 

It is interesting to observe that the amplitudes in (|5.29p are both given in a factorized form and they are related by 
the replacement d ^ d. 

We are now ready to evaluate the concurrence C in full order. The relevant matrix elements in (j2.4p are calculated, 
for a (spatially) long (i.e., an almost monochromatic) incident wave packet who ^ 1, impinging on a "small" target 
w ^ d, to be 

an = / d^^Kfeitil^lV^o)!' 
J An 

re-['''''^''"'''l/'"''|AA,j2|„,^_e''=-Xl' / d'fc |1 - e-''=°^"^(l - 6o)/feol', (5.33a) 



27r3;i4w2 

022= / d^fcKfciitl^lV'o)!' 

J An 

^2 



An 



-[d--(fco-d) _ e-^'^-XP / d^k |1 - e''=«'="*(l - 60)^0^ (5.33b) 



ai2= / d^k{kin\S[M{kiinS\^^o)* 
J An 



27r3;i4w2 

rf2feg.feofc-d[i _ g-.fcofc-d(i _ _ e*'^"*(l - CfeJ/feJ*. (5.33c) 



AO 
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We arrive at the following expressions for the concurrence 
2|ai2| 



C{p) 



ail + 022 



2|(afc„- 








6o)A.o][l-e'^°^"^(l-a-„)/feo]* 






Xr / d'fc|l-e-*'=°^"'(l-efco)/feor + |afco-e 


-iko-di 


d2fc|l-e^'=°^'*(l-ao)/fco|' 



(5.34) 



and the yield 

P = an + a22 



m 



-[d''^(ko-df]/Siu''\j^^ 



Ofco - e 



ika-di |2 

"fco I 

— ihc-d-' I 



d2fc|l-e-*'="^"*(l-a-o)/feo| 



J An 



(5.35) 



in full order. Observe that the dependence on the incident angle 0o remains through the quantity e 
this case with w ^ d, which is nearly absent at the lowest order [comment (i) in IIV A l] . On the other hand, the 
dependences of the concurrence C and of the yield P on the scattering angle Ob appear only through the quantity 
(see Appendix El 



C 



An 



^ikok d ^ ^ikgdcosBn 



An 



1 - -(kod)'^ sin^ sin2(A6l/2) - ikodcoseu sin2(A6l/2) 



(5.36) 



the absolute value of which is nothing but the concurrence at the lowest order (|4.15|) and one can further reduce C 
and P to obtain 



Cip) 



y 



where 



X ^ 2ia,, - e'^"-X)K„ - e"'"-%o)iC - (1 - Ck„)fk,, - (1 - 6o)*A; +C*\1~ ikof\fko\'] 

y ~ - e*'-xni + ii - ^kjVkoi' - c(i - (k„rn„ - c*(i - ajAoi 

+ |afc„ - e-^^-'^fefe^n^ ^ 1^ _ efcoH/fcoP - C(l - CfcJA,, - C*(l - aj*A;]- 



(5.37) 

(5.38a) 
(5.38b) 



In Fig. [6l the concurrence C and the probability P in full order are presented. Due to the resonant scattering of X 
between A and B, the concurrence C oscillates as a function of the scattering angle Ou [Fig. [Hlja)]. Such an oscillation 
is absent at the lowest order in gr (Fig. [2]). The resonant scattering can enhance the entanglement, if X is captured 
in appropriate directions [l3|. The overall characteristics (i)-(iv) of the concurrence C discussed in Sec. IIVI for the 
lowest-order estimation (Fig. [2]) are however more or less kept, apart from the oscillation around the lowest-order 
values. As shown in Fig. [7]^a), the effect of the resonant scattering saturates for mgr/h^d ^ 1. In this regime, the 
concurrence C and the probability P in (|5.37p are given in terms of 



y- 



'2iki)d 



'ikod 



Ah(-)-d 



ikod 



ikod 



-e 



iko-d 



kid^ 



ryikad 



iko-d 



ikod 



1 / *e 

- 2Im C* 



kld^ 



C-2 



ikod 



sin fend C* 



kod 



kld^ 



(5.39a) 



k^d 



^2ikod 



kld^ 



1 - 



iknd 



^ — ikn-d 



kld^ 



ikod 



kod 



(5.39b) 



and 



im-^grkl (1 + e^^'^'xi/k^d'^y 



(5.39c) 
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FIG. 6: (Color online) (a) Concurrence C and (b) yield P in the full order of perturbation with a large incident wave packet, 
w ^ d, shown in polar coordinates with radii C and P as functions of the scattering angle relative to the alignment d of the 
target qubits A and B. Parameters are: kod = 10, 9o = 90°, AS = 7r/15, and mgr/h^d — 1. The corresponding concurrence C 
in the Born approximation is also shown in (a) as a reference (dashed curve). 




FIG. 7: (Color online) (a) Concurrence C and (b) yield P in the full order of perturbation with a large incident wave packet, 
w ^ d, as functions of the scattering angle 9n and the coupling constant Qr- Parameters are: kod = 10, 6o = 90°, and 
Ae = 7v/15. 

Differently from the lowest-order result, the oscillation of the probability P (Young- like interference) is observed in 
Fig. El^b) [cf. Eq. (|4.14[) and Sec. lIV A^l for the lowest-order estimation]. This is because the multiple scattering with 
spin flips breaks the one-to-one correspondence between the path taken by X and the spin state of A and B after 
the scattering. Since the two paths are less distinguishable than at the lowest order, they interfere. As the coupling 
constant gr is increased, the probability P for the detection of X in the down state ||)x grows for mgr/h^d < 1, while 
it is suppressed for mgr/ti^d ^ 1 [Fig.[7jb)]. Sec (|5.39|) for the mgr/h^d » 1 case. The probability decays hke ~ g~'^ . 

As already mentioned above, the concurrence C and the probability P in full order depend on the incident angle 6*0 
even for a large incident wave packet w ^ d. See Figs. [5] and [5) The dependences are more prominent for a smaller 
incident wave number fco and a larger coupling constant gr- 
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FIG. 8: (Color online) Concurrence C in the full order of perturbation with a large incident wave packet, S> rf, as functions 
of the scattering angle 9t> and the incident angle 9o. Parameters are: (a) kod = 10, mgr/h^d = 1; (b) kod = 10, mgr/h^d — 10; 
(c) kod — TV, mgr/h^d = 1; (d) kod = tt, mgr/h^d = 10. AG = -it/15 for all cases. 




FIG. 9; (Color online) Probabilities P corresponding to the concurrences in Fig. |8] 
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VI. SUMMARY AND OUTLOOK 



111 this paper, a two-qubit system A+B, which is prepared initiaUy in a product state and has no mutual interaction, 
is shown to be made (maximally) entangled when another mediator qubit X is found spin-flipped in proper directions 
after scattering. The interaction between the mediator and each of the target qubits is described by a delta-shaped 
potential multiplied by the spin-exchange operator. The scattering (S) matrix elements are calculated both to lowest 
order and in full order in perturbation to estimate the concurrence C of the target system A-l-B and its yield P. 
The concurrence C, a measure of the entanglement, is evaluated as a function of such parameters as the momentum 
hko, width of the wave packet w of the incident qubit, the scattering (detecting) angle 0d and the detector resolution 
A0. The condition to obtain higher entanglement (concurrence) is derived as a function of these parameters and is 
interpreted in the context of the (in)distinguishability of alternative paths for quantum particle. Even though it is 
based on the lowest-order expression of the concurrence ()4.13p or (j4.17p . the condition seems to be 'globally' valid in 
full order in perturbation, except 'local' oscillations that are due to the presence of multiple scatterings. In particular, 
we may conclude that a highest (or maximal) value of concurrence C can be expected in the target system whenever 
the mediator qubit X is endowed with a short wave length fcgw ^ 1 and a spatially long wave packet w ^ d (compared 
with the size of the target d) and is injected and detected by a detector with a higher resolution A9 <C 1 on the line 
extending the alignment d of the target qubits A and B, that is, = 6*0 = 0. 

The result, i.e., the maximal entanglement (i.e., the largest concurrence C = 1) can be obtained in the target 
qubits A-l-B whenever the qubit X is injected and captured in the direction d of A-l-B, is considered to be remarkable 
since we do not need to adjust any parameters to obtain the maximally entangled state, contrary to the cases with 
controllable interaction times, where the acquisition of maximal entanglement is conditioned to the proper adjustment 
of the interaction strength, and to the cases of one-dimensional scattering, where the resonant effects play a crucial 
role Q . Similar results have also been obtained in the two-dimensional case [l^ . One could think that in the direction 
parallel to d, the lowest-order contribution, where the path-indistinguishability becomes maximum, resulting in the 
maximum entanglement, could always overwhelm the higher-order ones reflecting the multiple scatterings. This does 
not happen in one-dimensional cases, where there are no other spatial degrees of freedom for the particle to escape, 
while in spatial dimensions greater than one, the particle could be scattered in other directions than d, reducing such 
resonant effects on the concurrence and keeping the maximal concurrence. 
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APPENDIX A: ANALYTICAL ESTIMATION OF THE CONCURRENCE C AT THE LOWEST ORDER 

The scattering-angle dependence of the concurrence C and the yield P appears essentially through the quantity 



Acknowledgments 




(Al) 



the absolute value of which is nothing but the concurrence at the lowest order [see (|4.15p ] and can be evaluated as 
follows. First, the limitation on the polar angle variable < 9 < A6, where 9 is measured from the direction of the 
detector (that is, the scattering direction) fc^ and AQ = 47r sin^(A6'/2), is cast into an integration form 
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Now the two-dimensional angle integrations over k arc trivially done to yield 

dX- ^ 

oo {\ — ie)\kod + Xk-ol 

kod y„oo - cos Ob - ie) VC^Tsh?^ 

where we have changed the integration variable to ^ = X/{kod) + cos^d- Notice that this expression is exact without 
any approximation. 

Since either exponent of the integrand cannot be well approximated as a quadratic function, in particular for the 
case of small A^, in which case the exponent varies quite slowly and higher-order terms could not be neglected, the 
stationary-phase approximation does not work. We instead evaluate the integral (|A3[) in the complex ^ plane. First 
observe that each term of the integrand has a simple pole at ^ = cos 9d + ie in the upper-half plane and two branch 
points at ^ = ±i sin^D- Draw the two cuts from these branch points, parallel to the real ^ axis, to plus-infinity. Notice 
that the quantity \/ S,^ + sin^ Ou = \/ {£, + i sin 6*0) — i sin 6'd) behaves like 



e + sin^ eu ^ -l^le'"^ , = Arg(e + i sinfe), (A4) 

at infinity |^| — > 00, since the argument tp is measured anti-clockwise from the positive real axis, < f < 2tt, and 
then that of ^ — i sin^D is expressed as — (27r — ip) at infinity, except at infinity on the real positive ^ axis, where the 
arguments are both zero. The choice of the arguments is consistent with the fact that the above quantity becomes 
positive on the entire real ^ axis, where the arguments of ^ -|- i sin^D and S, — i sin^D are ip and —ip and cancel each 
other, as they should be. The argument in the region surrounded by the two cuts in the right-half plane will be 
scrutinized separately later. Therefore the first term in (jA3|) can be evaluated by deforming the integration contour 
downward, for it vanishes at infinity in the lower-half ^ plane. It is given by the integration along a contour that 
starts from +00 — i sin^D and runs first to the branch point at ~i sin^D below the (lower) cut, then to the origin and 
finally to +00 on the real positive ^ axis. It becomes, apart from the factor —l/kod, 



f 

Jo 



> ^—ikod{x — e^^^ ) cos AO ^ — ik(jdy^x{x — 2i sin 6'd) /"Oo ^ — ikod{x — cos 60) cos Ad ^ht^d^J a;2-|-sin^ 

dx rr, , + / dx ■ 



X - e^^D - ie ^x{x - 2i sin 6'd) h cos 6'd - ie ^ x"^ + sin^ 6'd 



/O g— 2fcocZ(iy— cos^d) cos gifcocJ-^sin^ 9u — y'^ 
dy — (A5) 
sinSo z2;-cos6'D-ze v^sin^ 6'd - 



The second term can be evaluated by similarly deforming the contour but this time upward on the complex ^ plane. 
The expression reads, again apart form the same factor —l/k^d, as 



! ^—ikQd{x — e ^^^] cos AO ^^kod^y x{x+2i sin 9u) ^—ikod{x— cos Ou) cos ^—ikod^J x^+sin^ Ou 

dx ■ , — / dx ■ 



x-e ^f^o - ie ^x{x + 2i sin 6'd) Jq x- cos 0^ - ie ^ x"^ + sin^ 6Id 

^sin p — ikQd{iy— cos ^d) cos ^—ikody^sin^ Ou — y^ 

+ i dy ■-■'^^,>"'fco''fl(-r.nsflp) (A6) 

Jo W - cos 6Id - ie \/sin^ 6'd - 

Notice that as for the first terms in the above expressions, representing the contributions coming from integrations 
along the lower and upper cuts, the integration contours can further be "pushed" up (down) by sin (in other 
Riemannian sheets), since the argument of x{x — 2i sIu^d) [or -^Z x{x + 2i sin ^d)] remains negative (or positive) 
there. That is, we have 



oo ^ — ikod[x — e^^^) cos AO ^—ik^d^ x(x — 2i sin ^?d) 

dx 

X - e»«D - ie v^a;(x-2isin6'D) 



roo ^—ikod{x— cos Oj:,) cos A9 ^ — ikady/ x'^-^sin-^ 9u rO ^—ikQd{iy— cos 9u) cos A9 ^ — ik od-y^ sin-^ OY) — y'^ 

I dx , =r- +i dy ^ =- (A7) 

Jo X- cos 6'd - le ^a;2 4, gj^^ 6Id J- sin Ou i-V " cos 6Id - le ^gin^ _ y2 
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FIG. 10: (Color online) Integration contours. The original contour Co and the deformed ones C±. 
(evaluated along the contour C- in Fig. [T0|) and (along the another contour C+) 



dx 



^—ik^d[x — e ^^^) cos A9 ^ikf)d^x{x+2isin9j:,) 



— cos ) cos A6' ^ikQd-\J a;^H-sin^ 

x-cos6iD-ie Va;^ + sin^ 6*0 



a; — e — ie 



/•sin Sd g 

i / dy- 



-ikQd{iy~cos 61^) COS A9 ^ik^dy^ sin^ Ou — y^ 



ly 



+ 2TTie''""^0{coseu). 
We thus obtain the following (still exact) expression for C Ail, 

CAn = -2ni——e{coseu) + 2'Ki , , 6>(-cos6'd) 



- cos Ob - ie a/sih^ 9d - y^ 



kod 

i , e 

- T—j dy 



kod 

— ikod^iy — cos Oj:) J COS L^u g 

— cos 9t) — ie 



cos Af? ^ikodyf sin^ Oiy—y'^ _|_ g— */cnf^'\/siii^ ^D^y^ 



V sin^ - 2/^ 



The last integration over y, which is rewritten as an integral over angle through y — sin^D sine 



/'37r/2 gi/corf sin 0d (cos 0— i cos A^? sin <; 



i sin sin (j) — cos 0d 



pi/corf cos COS A0 



can 



1 be cast into an integration along a unit circle ^ = e"^ enclosing the origin 
or equivalently that over z = 1/Q running clockwise, 

9 r ^ikodsinGii cos'^ {A9 /2)z ^ikodsin Oj:, sin^ {A6 /2) / z 

^ • ■ ifcodcos 6Id cos AS 

\z + cot(6lD/2)] [z - tan(6'D/2)] 



(A8) 



anti-clockwise on the complex 



kod sin 9u 



I 

dz ■ 



KoC 

and evaluated as the (infinite) sum of the residues within the unit circle. 



(A9) 

(AlO) 
plane, 

(All) 



gifcod g — ifco 

2TTi— — -0(cos 9b) — 2'Ki— — - 
kod kod 



— ikod 



(-cos 6*0) 



^ [ikodsm9Dsm^{A9/2)]"'+'^ 1 
fcodsin6'D ^ (n + 1)! 

n— ^ ^ 



gifcoc/ sin cos^(A^/2)^ 

n! rfz" [z + cot(6'D /2)] [z - tan(6'D /2)] 



ikod cos cos 
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= 27ri— — -6'(cos ^d) — 27ri— — —O^— cos ^d) 
kod kod 

- 7^ E E , , ,„„ |[2iMcos^(gD/2) sin2(A0/2)]"+i[2iWsin2(eD/2) cos2(A0/2)]^ 

- [- 2ifcodsin2(6'D/2)sin2(A6l/2)]"+i[-2ifco(i 0082(6*0/2) cos2(A6'/2)]^}e'''^«''^°'^''°™^^^ 

(A12) 

Therefore, we arrive at 

- oo n ^ 

C An = --^ y V- — ([2iMcos2(0D/2) sin2(A0/2)]"+i[2iMsin2(0D/2) cos2(A0/2)]^ 

knd ^ — ' ^ — ' n + ^. 

^ n=0 £=0 ^ ' 

- [-2ifcorfsin2(6lD/2)sin2(Ae'/2)]"+i[-2i/fcodcos2(6lD/2)cos2(A6l/2)]^}e''*=°''™"®°'=°"'^^ (A13) 

This is an exact expression, though the double summation seems to be difRcult to perform. On the other hand, it is 
suited for expansion for smaU A0. Indeed, it aUows us to obtain an estimation of C Afi, up to, say order (A^)^, 



C AO - [2ifcodsin2(A6'/2) - (fcod)^ sin^(A6'/4)(2cos6'D + ifco^^sin^ 6*0)] e*"i'^°'^«D 



kod 

1 



1 - -(fcod)^sin^ Ob s\tl\M / 2) - ik^d cos s\tl\M/2) 



(A14) 



which is (|5.36p and also yields an approximate expression for the concurrence C to lowest order (|4.17p . 

C = |C| - 1 - ^(fcod)^ sin^ 6*0 sin2(A6'/2). (A15) 

APPENDIX B: PROJECTION OPERATORS FOR THREE-QUBIT SYSTEM 

A quantum system composed of three spin-1/2 particles, say X, A, and B, is decomposed into a spin-3/2 and two 
spin-1/2 sectors, 1/2 (X) 1/2 (g) 1/2 = 3/2 © 1/2 © 1/2. The spin-3/2 states are totally symmetric with respect to the 
exchange among three particles, while the spin-1/2 states have mixed symmetries and one of them is anti-symmetric 
(singlet state) and the other symmetric (triplet state) under the exchange between, say A and B. The projection 
operators Vsij), Qs{l), and Ti.± introduced in Sec. IV Bl are related with these spin decompositions. 

First of all, since the square of the total spin is expressed as 

Sl^^ = ^ + <t(A) + ^(B)) ' = ^ (9 + 2^(^) • (<t(A) + ^(B)) + 2^(A) • <t(b)) , (Bl) 

the operator 7?.+ extracts the total spin 3/2 states, while TZ- is a projection upon the four-dimensional subspacc with 
the total spin 1/2. That is, 

3 1 

TZ+ = 1 <^ total spin -, TZ+ = total spin - (B2) 

or 

1 3 

7e_ = 1 <^ total spin -, 7e_ = <^ total spin -. (B3) 

Since Qs{j) is orthogonal to TZ+, while the subspace projected by the latter is included by that by Vsij) [see 
Eq. (|5.2ip ]. the projection operator Vsij) extracts all spin-3/2 states and a part of spin-1/2 states. Indeed, it is 
not difRcult to show that the latter sector with the total spin 1/2 (belonging to TZ- = 1 sector) are spanned by the 
following eigenvectors 

3 



|V^+) = yT^|0+)xAB-^s^^|T)x|O,O)AB, IV--) = v/i^I'/'-)xab-*s^^U)x|0,0)ab, p - - - 37, (B4) 
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where states \(j)±)xAB are particular states belonging to the AB-triplet subspace 

I0+)XAB = y^|T)x|l,0)AB - y|u)x|TT)AB = (|T)a|0, 0)xB + |T)b|0, 0)xa) , (B5) 

|'/'-)xAB - y||T>xUi)AB - \/|u)x|1,0)ab = -^(U)a|0,0)xb + U)b|0,0)xa), (B6) 

and |0, 0)xA the singlet XA-state, and so on. That is, wighted superpositions of the AB-triplet and AB-singlet states, 
characterized by the parameter p (or 7), as well as the total spin 3/2 states, are extracted by the projection operator 
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